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ABSTRACT 

Classical novae are stellar explosions occurring in binary systems, consisting of a white dwarf and a 
main sequence companion. Thermonuclear runaways on the surface of massive white dwarfs, consist- 
ing of oxygen and neon, are believed to reach peak temperatures of several hundred million kelvin. 
These temperatures are strongly correlated with the underlying white dwarf mass. The observational 
counterparts of such models are likely associated with outbursts that show strong spectral lines of neon 
in their shells (neon novae). The goals of this work are to investigate how useful elemental abundances 
are for constraining the peak temperatures achieved during these outbursts and determine how robust 
"nova thermometers" are with respect to uncertain nuclear physics input. We present updated ob- 
served abundances in neon novae and perform a series of hydrodynamic simulations for several white 
dwarf masses. We find that the most useful thermometers, N/0, N/Al, 0/S, S/Al, 0/Na, Na/Al, 
0/P, and P/Al, are those with the steepest monotonic dependence on peak temperature. The sensitiv- 
ity of these thermometers to thermonuclear reaction rate variations is explored using post-processing 
nucleosynthesis simulations. The ratios N/0, N/Al, 0/Na, and Na/Al are robust, meaning they are 
minimally affected by uncertain rates. However, their dependence on peak temperature is relatively 
weak. The ratios 0/S, S/Al, 0/P, and P/Al reveal strong dependences on temperature and the poorly 
known ^^P(p,7)^^S rate. We compare our model predictions to neon nova observations and obtain 
the following estimates for the underlying white dwarf masses: 1.34-1.35 (V838 Her), 1.18-1.21 
Mo (V382 Vel), <1.3 M© (V693 CrA), <1.2 M© (LMC 1990#1), and <1.2 M© (QU Vul). 
Keywords: Stars: abundances, novae, cataclysmic variables, nuclear reactions, nucleosynthesis, abun- 
dances 
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1. INTRODUCTION 

Classical novae are stellar explosions occurring in bi- 
nary systems co nsisting of a white d warf and a close 
companion star (jBode fc EvansI l2QQ8f ). When the stel- 
lar companion expands beyond its Roche Lobe, a mass 
transfer episode ensues. This matter possesses angular 
momentum due to the system's rotation and forms an ac- 
cretion disk, falling onto the surface of the white dwarf. 
At some point, matter from the underlying white dwarf 
is mixed with the accreted material, and this mixture 
is heated and compressed until a thermonuclear runaway 
takes place, ejecting matter into the interstellar medium. 
The ejected material consists of a mixture of white 
dwarf and accreted materials that h as been processed 
by explosive hydrogen burning (Starrfield. Iliadis. fc Hixl 
l2QMlShaxa et al.lllmllYaron et al.ll2QQ5L and references 
therein). Spectroscopic studies have identified two dis- 
tinct white dwarf compositions, carbon-oxygen (CO) and 
oxygen- neon (ONe). Nova shells rich in CNO material 
point to an underlying CO white dwarf, which represents 
the evolutionary fate of a low-mass star after core helium 



burning. On the other hand, element enrichments^ in the 
range of Ne to Ar have been attributed to underlying, 
massive ONe white dwarfs {Mwd ~ 1-1 M©), represent- 
ing the evolutionary fate of stars ( M ^ 9-11 M^;^) after 
completion of core carbon burning (jAlthaus et al.|[2QTQl 
and references therein). The latter explosions, which are 
often referred to as neon n ovae, tend to be much more 
energetic than CO novae ( Starrfield. Sparks, fc TruranI 
1986) and are the main interest of the present work. 

The study of classical novae is of significant in- 
terest for several reasons. Spectroscopic studies, af- 
ter proper interpretation, constrain stellar evolution- 
ary models by revealing the composition of the un- 
derlying white dwarf. Observed elemental abundances 
also provide a record of the thermonuclear runaway, 
including peak temperatures and expansion timescales, 
and have been used to co nstrain stellar explosion mod- 
els ([Starrfield eFaH l2QQQ[ ). Classical novae are note- 
worthy contributors to the chemical evolution of the 
Galaxy because they are major sources of ^^C, -'^^N , 
and 1^0 (|Kovetz fc Prialnikl[T997l : [jose fc Hernanzlll998l : 

^ Neon novae are mainly identifie d by the e arly appearance of 
the [Ne II] 12.8 ^im emission fine (,Gehrz et al.iri998l 
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Jose, Hernanz fc IliadisI I2QQ6I : iStarrfield, Iliadis, fc Hixl 
2QQ8I ) . Significant efforts have been undertaken to detect 
7-rays from novae as well. Although no nuclear 7-rays 
from classical novae have b een observed by any satellite 
observatory (Hernanz 2008), it has been proposed that 
7-rays should be present in these novae due to the decay 
of radioactive ^^F. 

Although theoretical models of classical novae are 
successful in reproducing the overall characteristics of 
the observed outbursts, several key issues remain un- 
explained and call for improved observations, simula- 
tions, and nuclear laboratory measurements. First, the 
masses of ONe white dwarfs are poorly constrained by 
observations. Models of intermediate-mass stars pre- 
dict val ues in the ran^e of ^ 1.1-1.4 (see discus- 
sion in iDohertv et al.l |2QTq1 ). Clearly, accurate white 
dwarf mass determinations are highly desirable. Sec- 
ond, hydrodynamic simulations underpredict the ob- 
served eject a masses by an ord er of magnitude or more 
(|Starrfield. Iliadis. fc Hi^l2QQ8f ). Third, although it is 
sometimes claimed in the literature that theoretical mod- 
els more or less reproduce the observed elemental abun- 
dance patterns of classical novae, we will show that such 
claims are premature regarding neon novae. 

The goal of the present work is to investigate if and how 
simulated elemental abundances in classical nova shells 
can strongly constrain important explosion parameters, 
most notably the peak temperature. In particular, we 
seek a robust relationship between model peak tempera- 
ture and simulated elemental abundance ratios. Such a 
classical nova thermometer^ if it exists, should be useful 
for improving models of thermonuclear runaways. Reli- 
able determination of the peak temperature during the 
outburst will also greatly improve estimates of the mass 
of the underlying white dwarf. 

In this work, we focus on exploring the peak tempera- 
ture (or, equivalently, white dwarf mass) parameter space 
of classical nova models. Clearly, the degree of mixing 
between white dwarf matter and the accumulated mat- 
ter before the explosion, or m ore general ly, the mixing 
mechanism (see discussion in Casanova et al. 2011), sig- 
nificantly impacts the ejected abundances. However, in 
this work we will assume a mixing fraction of 50% and 
will leave a study of varying this parameter to future 
work. 

From a nuclear physics point of view, the differences 
of the nucleosynthesis pattern in the A < 20 and A> 20 
mass regions are readily apparent. In the former case, the 
hot CNO cycles give rise to a redistribution of material 
in C, N, and O and, in favorable cases, the peak temper- 
ature is likely reflected in the observed abundance ratios. 
In the latter case, the nucleosynthesis is characterized by 
a continuous flow from Ne up, perhaps infl uenced by re - 
action cycles at ^^Na, ^^Al, ^ip^ and ^^Cl miadisll2007D . 

Thermonuclear reaction rates are highly sensitive to 
the stellar temperature and, consequently, small temper- 
ature changes have a large impact on elemental abun- 
dance ratios. This argument represents our starting 
point for studying potential thermometers for classical 
novae. For example, ^^Ne(p,7)^^Na is among the slow- 
est nuclear reactions during the thermonuclear runaway. 
The main reason is the low reaction Q- value, giving rise 
to a very small number of resonances at low center-of- 
mass energies. In fact, the large Ne abundance observed 



in the ejecta of neon novae is mainly caused by the low 
proton capture rate on ^^Ne^. 

Observed and predicted classical nova abundances have 
been comp ared in many work s over the years (see, for 
example, Starrfield et al. 1998; Kovetz fc PrialniWigal 
Jose & Hernanz 1998, and references therein). Models of 
classical novae for a ran ge of ONe whi t e dwa rf masses 
have been presented in iPolitano et al.l ()1995f ). but no 
attempt was made in that work to predict quantita- 
tively the white dwarf mass of a given observed neon 
nova. It should also be mentioned that these models 
did not result in mass ejection, so their results must be 
treated wit h care. Nova model sequences have als o been 
evolved by IWanajo. Hashimoto, fc Nomotol ()1999f ) , with 
the aim of constraining ONe white dwarf masses of ob- 
served neon novae. However, this work does not present 
self-consistent models and must also be handled care- 
fully. As will be seen later, some of their predictions 
differ vastly from our results. Clearly, the nova models 
presented in the latter two works are based on thermonu- 
clear reaction rates that are about 20 years old and thus 
are completely outdated. 

Furthermore, it must be emphasized that none of the 
studies to date explored the crucial effects of thermonu- 
clear reaction rate variations on the predicted abun- 
dance ratios. A key aspect of our work is the use of 
a newly a vailable nuclear physics library, called STAR- 
LIB (Ilia dis et al.l [2011). It incorporate s a recent eval- 
uation of thermon uclear reaction rates (Longland e t al.l 
I2OIOI : llliadis et all L2010a.b.c), obtained from a Monte 
Carlo sampling of experimental nuclear physics observ- 
ables, and provides for the first time statistically mean- 
ingful reaction rates and associated uncertainties. 

In Sec. 2, we present updated observed abundances 
in neon nova shells. It will become clear that it is ad- 
vantageous to focus the discussion on elemental abun- 
dance ratios rather than the abundance of one element. 
Our methods and results are presented in Sees. 3 and 
4, respectively. First, we computed four new hydrody- 
namic models of classical novae, with peak temperatures 
in the range of 228-313 MK. The subset of elemental 
abundance ratios revealing the steepest dependence on 
peak temperature, that is, with the promise of poten- 
tial use as a thermometer, was considered for further 
inspection. Then, suitable post-processing calculations 
were performed, demonstrating that the final predicted 
elemental abundances approximately agree with those 
from the self-consistent hydrodynamic model. In a fi- 
nal step, many post-processing network calculations were 
performed by varying one reaction rate at a time within 
its current uncertainty. The impact of this sensitivity 
study on potential classical nova thermometers was then 
investigated. Our results are compared to observation in 
Sec. 5. A summary and conclusions are given in Sec. 6. 

2. OBSERVED ABUNDANCES 

Before discussing our new simulations, we will 
summarize the abundances observed in neon nova 
shells. Compilations of classical nova abundances 
(expressed as mass fr a ctions ) can be found, for ex - 
ample, in ICehrz et al.l (jl998f ): iStarrfield et~al1 (jl998f ): 

^ The operation of the ^^Na(p,Q;)^°Ne reaction also contributes 
to the survival of ^°Ne. 
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IWanajo, Hashimoto, &: Nomotol (| 19991 ). First, these 
compilations reveal that in several cases, even for iden- 
tical events, different authors report rather different el- 
emental abundances, with differences amounting up to 
an order of magnitude for some elements. Second, no 
uncertainties are provided with the listed abundances, 
making it difficult to assess the significance of a given 
element observation. Third, some of the compilations 
provide only the mass fraction sum for elements between 
Ne and Fe. As will become apparent later, when mod- 
ehng neon novae the individual observed abundances of 
intermediate-mass elements (Ne, Mg, Al, Si, and S) are 
extremely important. Fourth, for several novae listed 
in the above compilations, the mass fraction sum differs 
significantly from unity, reflecting errors in the reported 
values. Since all of these reasons could contribute to the 
mistaken impression that reported nova abundances are 
highly uncertain and scatter widely, we felt compelled to 
present updated abundances. 

It is difficult to determine reliable abundances from ob- 
served nova ejecta. For example, the derived abundances 
are model-dependent since most analyses assume spher- 
ical shells although there is evidence that t he eje cta are 
not spherically symmetric (jCasanova et aL|[2Qll[ ). Fur- 
thermore, the ejecta are potentially chemically inhomo- 
geneous, adding substantially to the complexity of the 
problem. Also, the filling fraction (i.e., the fraction of 
the shell volume occupied by gas) is poorly constrained. 
Finally, the abundance analysis must account for the 
(sometimes substantial) fraction of unobserved ionization 
states. These, an d other, difficu l ties h ave been discussed 
in more detail bv lJose fc Shore! ()2QQ8f ). 

In the past different authors have used different pro- 
cedures in the analysis of optical, infrared, and ultra- 
violet spectra. However, over the past 15 years many 
nova ejecta have been analyzed using consistent tech- 
niques, and the situation has improved significantly. In 
particular, we list all known neon novae in Table [TJ with 
the exception of VI 370 Aql and U Sco^. For the eight 
neon novae listed, the abundances were analyzed using 
the photoionization code CLOUDY (Ferland 2003) cou- 
pled to the least squares minimization code MINUIT, 
adopted from the CERN library. In this manner, the 
intensities of all spectral lines are determined together 
via successive iterations during the minimization, while 
in addition uncertainties of the line intensities can be 
quantified in a meaningful way. 

Nova abundances are usually presented in the litera- 
ture as "relative to hydrogen relative to solar", i.e., the 
quantity {N,i/Nh) 

where N^i refers to the number abundance. However, the 
"solar" abundances have undergone signif icant revision 
over the past two decades (see discussion in Jose fc Shore! 
!2008f ) , so observations were normalized to the solar abun- 
dances of Lodders, Palme, & Gail (2009). In order to 
compare the observations to our simulations (discussed 

^ Abundance an alyses of V137Q Aql are presented i n 
ISniiders et a ll fT987l) and lA ndrea. Drechsel. Starrfleldl {1991), 
but they were not performed using iterative least squares minimiza- 
tion techniques and hence the associated uncertainties are difficult 
to assess. For U Sco, o nly a few abundances, such as O and Ne, 
are discussed in .Mason! (|2011l ). 



below) in a meaningful way, we first calculated, from 
the values of ^ reported in the literature, the quantity 
Nei/Nn and, subsequently, converted the results to the 
elemental mass fraction ratio, Xei/Xu- In some in- 
stances, the abundances presented in the literature had 
asymmetric uncertainties. For the sake of simplicity, the 
uncertainties in these cases were symmetrized by adopt- 
ing the central value between the upper and lower bounds 
as the recommended value. The resulting mass fraction 
ratios for elements between He and Fe are listed in Ta- 
ble [H The hydrogen mass fraction is then obtained from 

Xn = [l + ^ + ^ + ... + ^r (2) 

and is listed in the last row. 

The individual mass fractions can be easily calculated 
from the values of X^i/Xh and Xh- Note, however, that 
the values of X^i derived in this manner are sensitive to 
the abundance fraction missed in the spectral line anal- 
ysis, such as missing elements or unaccounted for ion- 
ization states. The ratio X^i/Xh is expected to be less 
susceptible to systematic errors than Xei and, therefore, 
we focus in the following sections on mass fraction ratios 
rather than mass fractions. 

Inspection of Tabled] reveals large overabundances rel- 
ative to the solar Xqi/Xh values, listed in the last col- 
umn. For some neon novae, N, O, Ne, Mg, Al, Si, S^ 
and Ar, are enriched by factors of 160, 50, 260, 40, 260, 
40, and 40, respectively. In the following sections, we 
will discuss the abundances predicted by models of neon 
novae. 

3. METHOD 

Our nuclear reaction network followed the evolution 
of 117 nuclides in the range from H to Ti through 635 
nuclear processes, including weak interactions, reactions 
of type (p,7), (p,Q^), (<^,7), and their reverse reactions. 
The rates of these interactions are a dopted from a next - 
generation library, called STARLIB (jlliadis et al.!!20lTh . 
A detailed account of this library will be given in a 
forthcoming paper. In brief, STARLIB is a tabular li- 
brary, listing the rates and corresponding uncertainties 
on a temperature grid between 1 MK and 10 GK for 
each nuclear interaction. For 62 of the nuclear reactions 
in the A = 14 — 40 mass range contained in STAR- 
LIB (among them most of the reactions important for 
nova nucleosynthesis), experimental rates are based on 
the 2010 Monte Carlo evaluation of thermonuclear reac- 
tion rates (jLondand et al.!!2010! : flliadis et al.!!2010a!lb! O). 
The Monte Carlo reaction rates and associated uncer- 
tainties are obtained from sampling randomly over un- 
certainties of (input) nuclear physics observables and are 

^ As will become apparent later, observations of sulfur are par- 
ticularly important. In addition to the novae listed in Table ^ 
sulfur has also been observed in the optical s pectra of QV Vul, 
V2214 Oph and V443 Set ( Andrea, Drechsel Sz Starrfield 1994). 
For these three objects, which are most likely associated with 
outbursts involving CO rather than ONe white dwarfs, the de- 
rived su lfur abundances are near the solar value. Furthermore, 
lAndrea. Drechsel, Starrfi eld (1994) report a relatively high sul- 
fur abundance for QU Vul. However, the optical spectrum shown 
in their Figure 1 reveals only a weak sulfur [S II] line. Their derived 
high abundan ce value is pe rhaps biased by their analysis method 
used. Finally, rSniiders (1987) report a high sulfur abundance 

for the neon nova V1370 Aql. Unfortunately, no uncertainty is pro- 
vided and the reliability of their reported value is difficult to assess. 
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Table 1 

Summary of Observed Abundances (in Mass Fractions) for Neon Novae from UV, optical, and IR Spectroscopy 





LMC 1990#li 


V4160 Sgr^ 


V838 Her2 


V382 Vel^ 


QU Vul^ 


V693 CrA^ 


V1974 Cyg^ 


V1065 Cen'^ 


Solar^ 


^He/^H 


4.8(8)E-01 


7.1(4)E-01 


5.6(4)E-01 


4.0(4)E-01 


4.6(3)E-01 


5.4(22)E-01 


4.8(8)E-01 


5.4(10)E-01 


3.85E-01 




3.7(15)E-02 


1.43(7)E-02 


2.28(23)E-02 


2.6(13)E-03 


9.5(59)E-04 


1.06(44)E-02 


3.1(9)E-03 




3.31E-03 




1.48(42)E-01 


1.27(8)E-01 


3.29(47)E-02 


2.28(54)E-02 


1.61(10)E-02 


1.84(67)E-01 


6.0(15)E-02 


1.40(33)E-01 


1.14E-03 




2.4(10)E-01 


1.35(9)E-01 


1.42(38)E-02 


4.13(38)E-02 


3.2(14)E-02 


1.63(66)E-01 


1.55(85)E-01 


4.7(15)E-01 


9.65E-03 




1.6(10)E-01 


1.38(5)E-01 


1.22(5)E-01 


4.0(7)E-02 


5.1(4)E-02 


6.7(34)E-01 


9.7(40)E-02 


5.34(98)E-01 


2.54E-03 




1.37(71)E-02 


?^8.4E-03 


1.2(7)E-03 


2.45(14)E-03 


1.02(49)E-02 


9(7)E-03 


4.3(28)E-03 


4.4(13)E-02 


9.55E-04 




2.3(ll)E-02 




1.8(13)E-03 


1.63(16)E-03 


4.1(ll)E-03 


5.0(46)E-03 


>7.8E-05 




8.74E-05 




4.8(39)E-02 


1.09(6)E-02 


7(2)E-03 


5(3)E-04 


2.4(18)E-03 


2.4(18)E-02 






1.08E-03 








1.48(15)E-02 










2.3(13)E-02 


5.17E-04 












4.0(3)E-05 






4.6(17)E-03 


1.29E-04 






2.4(8)E-03 


2.35(63)E-03 




9.53(54)E-04 




8.8(72)E-03 


1.16(40)E-02 


1.81E-03 




4.7(9)E-01 


4.65(37)E-01 


5.63(36)E-01 


6.6(4)E-01 


6.3(3)E-01 


3.8(14)E-01 


5.5(8)E-01 


3.6(10)E-01 


7.11E-01 



Note. — All abundances are given here in terms of mass fraction ratios, X^i/Xh (or mass fraction for hydrogen; see last row), by converting the "number abundances 
relative to hydrogen relative to solar" from the original literature (references provided below). The abundance uncertainties are given in parentheses. 

^ FromlVanlandingham et al.| (Il999f) . 
^ From[Schwarz et al. ( 2007). 
^ From|Shore et al. ( 2003^). 
^ From[Schwarz (2002). 

^ From^anlandingham, Starrfield . Shore! ^199% . 

^ From[yanlandingham et al. (200^; solar abundances assumed in their analysis are not listed; their adopted values were: log (A/'g^/A/'/f )0 = — 1.0 (He), —3.45 (C), —4.03 

(N), -3.13 (O), -3.93 (Ne), -4.42 (Mg), -5.53 (Al), -4.45 (Si), -4.79 (S), -4.49 fFe) d Vanlandinghamll2Qr2i) . 

^ From fHelTon et al.l (12^101) . 

^ From lLodders. Palme, Gaill (l2QQ9h . 

9 Calculated from: Xjj = [1 + Xhc/Xh -^Xc/Xh+ ••• +XFe/XH]~^. 
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Table 2 

Properties of Evolutionary Nova Models 



Property 




Model 




J115 


J125 


J 130 


J 135 




1.15 


1.25 


1.30 


1.35 


Tpeak (GK) 


0.228 


0.248 


0.265 


0.313 


Me," (lO-^Mo) 


2.46 


1.89 


1.17 


0.455 


R^'d (km) 


4326 


3788 


3297 


2255 



derived from the 0.16th, 0.50th, and 0.84th quantiles of 
the cumulative distribution of the (output) reaction rates 
(for a coverage probabihty of 68%). The information 
on rigorously defined reaction rate uncertainties was not 
available previously and opens interesting possibilities for 
nucleosynthesis studies. 

We generated a se ries of new hydrodynam ic models us- 
ms; the code SHIVA (Ij ose fc Hernanz|[l998f ) with reaction 
rates adopted from STARLIB. The different models were 
generated for a range of white dwarf masses (1.15 M© 
- 1.35 M0). Information on derived parameters of our 
models, specifically, peak temperature (Tpea/e), ejected 
mass (Mgj), and initial white dwarf radius (i?^*^), is 
given in Table O Each model included 45 envelope zones. 
Test calculations performed with 500 zones provided 
essentially the same results. The peak temperatures 
achieved in the hottest zone in our simulations, rang- 
ing from 228 MK to 313 MK depending on the model, 
can be regarded as typical for models of thermonuclear 
runaways involving ONe white dwarfs (for m odels that 
achieve higher peak temperatures, see Starrf ield et al.l 
[2009). We assumed a mixing fraction of 50% between 
accreted matter of solar composition and white dwarf 
matter prior t o the outburst. The sola r abun dances are 
adopted from iLodders. Palme, fc Gaill (|2009l ). while the 
latter abundances are taken from the evolution of a 10 
Mq star from the main sequence t o the end of core car- 
bon burning (Ritossa, Garcia-Berr o. fc Iben| [T996). The 
initial envelope composition for our nova simulations is 
presented in Table [3l The initial luminosity and mass 
accretion rate for all models amounted to Lini = 10~^ 
Lq and Mace = 2 x 10~^^ MQ/yr, respectively. 

Final isotopic abundances, for matter that reached and 
exceeded escape velocity (i.e., the fraction of the envelope 
effectively ejected), were determined 1 hour after peak 
temperature was achieved. By that time, short-lived 
parent nuclei had decayed to their stable daughters, and 
the isotopic abundances of a given element were summed 
for each model. Those elemental abundances revealing a 
strong, either monotonically increasing or decreasing, de- 
pendence on peak temperature were then considered for 
further study. 

Investigating the sensitivity of any useful nova ther- 
mometer to nuclear reaction rate variations implies 
that a large number of simulations need to be per- 
formed in order to exhaust the nuclear physics param- 
eter space. Clearly, full hydrodynamic model simula- 
tions are at present far too time-consuming for this 
purpose. Instead, we performed post-processing re- 
action network calculations, using time- temperature- 
density profiles adopted from the hydrodynamic models 
described above. Such calculations have been used suc- 



Table 3 

Initial Envelope 
Composition (Mass 
Fractions) of Present Nova 
Simulations 



Isotope 


Mass Fraction 


iH 


3.56E-01 


3 He 


4.23E-05 


4He 


1.37E-01 


6Li 


3.44E-10 


^Li 


4.91E-09 


9 Be 


7.51E-11 


lOg 


5.05E-10 


iiB 


2.26E-09 


12c 


5.74E-03 


13c 


1.42E-05 


14n 


4.04E-04 


15n 


1.59E-06 


I60 


2.59E-01 




1.37E-06 


180 


7.72E-06 


19p 


2.08E-07 


20 Ne 


1.57E-01 


2iNe 


2.99E-03 


22Ne 


2.22E-03 


23Na 


3.22E-02 


24 Mg 


2.77E-02 


25 Mg 


7.94E-03 


26 Mg 


4.98E-03 


27A1 


5.43E-03 


28Si 


3.51E-04 


29Si 


1.85E-05 


30Si 


1.26E-05 


31p 


3.50E-06 


32g 


1.74E-04 


33 S 


1.42E-06 


34 g 


8.24E-06 


35C1 


1.87E-06 


37C1 


6.29E-07 


36 Ar 


3.84E-05 


38 Ar 


7.40E-06 


39k 


1.86E-06 


40 Ca 


3.18E-05 


Note. 


— Values are 



obtained assuming 50% 
mixing of solar accreted 
matter with white dwarf 
material (see text). 

cessfully in pre vious reaction rate sensitivity studies of 
classical novae (|Iliadis et al.ll2QQ2l ). However, previous 
investigations employed one-zone calculations, appropri- 
ate for the hottest zone close to the white dwarf surface 
only. For the purposes of the present work it was im- 
portant to perform post-processing calculations for all 
burning zones^ matching the final abundances of the hy- 
drodynamic models as closely as possible. We explored 
two extreme assumptions regarding the mixing of ma- 
terial between zones to this end. In the first case, we 
performed post-processing calculations for each burning 
zone separately, using the appropriate time-temperature- 
density profiles from the hydrodynamic study, and at the 
end of the calculation the elemental abundances were 
mass- averaged over all zones. This assump tion is sim- 
ilar to the procedure used by Smith et al. (2002) and 
will be referred to as the "no- mixing approximation" . In 
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Figure 1. (Color online) Comparison of final elemental abun- 
dances from a full hydrodynamic model to results from post- 
processing calculations. Abundance ratios are displayed versus 
atomic number, Z. The nova model assumes accretion onto a 1.25 
M0 ONe white dwarf (Table O. The three panels are obtained for 
three different mixing assumptions in the post-processing calcula- 
tions. Top Panel: Instant-mixing; Middle Panel: No-mixing; 
Bottom Panel: Geometric mean of abundances from the top and 
middle panels. See text for details. Elements of interest in this 
work are displayed in color. Note the good agreement between 
hydrodynamic and post-processing results in the bottom panel. 
Similar plots are obtained for all other nova models explored in 
this work. 



the second case, the post-processing calculation is per- 
formed by replacing each local thermonuclear rate by its 
mass- weighted average over the convective region, im- 
plying that the turnover time is faster than the nuclear 
burning ti me. Th is assumption has been made, for ex- 
ample, bv iPrialnI k (1986) in the framework of a hydro- 
dynamic model and will be referred to in the following 
as the "instant-mixing approximation". The actual be- 
havior probably lies somewhere between these prescrip- 
tions. Therefore, a third set of elemental abundances 
was derived from the post-processing calculations, sim- 
ply by computing the geometric mean of the abundances 
resulting from the no-mixing and instant-mixing approx- 
imations. Representative results are shown in Figure [H 
comparing the final elemental abundances from a full hy- 
drodynamic simulation with the three sets of final post- 
processing abundances described above. It is especially 
encouraging that the geometric mean abundance values 



seem to reproduce the stellar model elemental predic- 
tions rather well especially for the elements of interest in 
this work. 

More than 7,000 multi-zone network calculations for all 
four classical nova hydrodynamic models were performed 
by independently varying the rates of 214 reactions, that 
is, the most important nuclear reactions in our network. 
For each reaction, the rates were multiplied by factors of 
100, 10, 5, 2, 0.5, 0.2, 0.1, and 0.01 in successive network 
calculations. Final elemental abundances were adopted 
from the mixing prescription (instant-mixing, no-mixing, 
or geometric mean) that best approximated the results of 
a given hydrodynamic model using recommended reac- 
tion rates. Final elemental abundances are then linearly 
interpolated between the above variation factors to deter- 
mine the impact of the actual reaction rate uncertainties, 
which were adopted from the STARLIB library. 

4. RESULTS 

We start with the results from our hydrodynamic sim- 
ulations using the code SHIVA. Figure [2] displays the 
final elemental abundances, normalized to initial abun- 
dance adopted in the simulation, versus atomic number, 
Z, for all four nova models (see Table [2j). It is apparent 
that certain elements are overproduced (N, Si, P, and S), 
while others are depleted (O, Na, Mg). We are particu- 
larly interested in final elemental abundances that show 
a steep dependence on peak temperature. We find that 
N, Na, P, and S monotonically increase with peak tem- 
perature while O and Al decrease. Note that F, CI, and 
Ar also show strong trends, but their absolute predicted 
final abundances are very smah (i.e., 10"'^ to 10"^ by 
mass). Therefore, they were deemed not useful for our 
present purpose. Based on the hydrodynamic simula- 
tions, our first indication of useful thermometers, with a 
steep dependence on peak temperature, are the eight ele- 
ment ratios N/0, N/Al, 0/S, S/Al, 0/Na, Na/Al, 0/P, 
and P/Al. 

The abundance ratios for these element pairs are shown 
in Figure [3] versus peak temperature, where the solid lines 
indicate ratios for elements that have been observed in 
neon nova shells (N, O, Al, S), and dashed lines denote 
ratios involving elements (Na, P) for which reliable abun- 
dances have not been reported in shells of classical no- 
vae yet^. The elemental abundance ratios 0/S (decreas- 
ing) and S/Al (increasing) display the strongest varia- 
tion between peak temperatures of 228 MK and 313 MK, 
amounting to about three orders of magnitude. Follow- 
ing closely are the abundance ratios 0/P (decreasing) 
and P/Al (increasing), which vary by more than two or- 
ders of magnitude. The other ratios shown in Figure [3] 
(0/Na, N/Al, N/0, Na/Al) show variations by about 
one order of magnitude. The next step was to investi- 
gate how robust these element ratios are with regard to 
thermonuclear reaction rate variations. 

The sensitivity of the nova thermometers to the nuclear 
physics input is shown in Figure IH Elemental abundance 
(mass fraction) ratios are displayed versus peak temper- 
ature, for the four nova models explored here (see Ta- 
ble [2]). Shown on the left-hand side are abundance ratios 

^ Infrared coronal lines from phosphorus, [P VIII] and [ P VII] , in 
the ej ect a of V1974 Cyg have been reported by Wagner DePovl 
but no elemental abundances have been reported. 
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Figure 2. (Color online) Final elemental abundances from all four 
hydrodynamic models explored in this work, normalized to initial 
abundances, versus atomic number. Some elements are overpro- 
duced (N, Si, P, S), while others are underproduced (O, Na, Mg). 
Of interest here is the monotonic dependence of certain elements 
on the white dwarf mass or, equivalently, peak temperature, re- 
sulting in a color sequence of red, blue, green, black, or vice versa. 




Peak temperature (GK) 

Figure 3. (Color online) Ratios of 8 elemental abundances (mass 
fractions), derived from hydrodynamic nova models, that show a 
steep and monotonic dependence on peak temperature, i.e., by at 
least an order of magnitude. These ratios are prime candidates for 
nova thermometers. 

involving elements that have been observed in neon nova 
shells, while the right-hand side displays abundance ra- 
tios involving at least one element that has not yet been 
observed. The solid red lines correspond to final abun- 
dance ratios obtained from the full hydrodynamic mod- 



els, while the solid black lines show the results of the 
mixing prescription that best reproduced the hydrody- 
namic model results using recommended reaction rates. 
Note that the red and black solid lines generally differ 
in magnitude by less than a factor of 2, providing sup- 
port for our conjecture that the present post-processing 
simulations reliably approximate the results from the full 
hydrodynamic models for the elements of interest in this 
work. 

The dashed lines in each panel of Figure 2] display the 
change of predicted elemental abundance (mass fraction) 
ratios caused by varying individual reaction rates within 
their unc ertainties (defined by a coverage probabil ity of 
68%, see iLondand et al.llMol : llliadis et al]l2QlQaD . In- 
terestingly, four of the thermometers (N/0, N/Al, 0/Na, 
and Na/Al) reveal rather robust abundance ratios, with 
an uncertainty of less than 30%. Consequently, these ra- 
tios represent currently useful nova thermometers, with- 
out the immediate need for improved nuclear laboratory 
measurements. However, only the first two (N/0 and 
N / Al) involve elements that have so far been observed in 
neon nova shells. The last two involve sodium, which has 
been se arched for in infrared spectra but not yet detected 
(IShorelfMl ). 

The use of the other four predicted elemental ratios 
(0/S, S/Al, 0/P, and P/Al) as nova thermometers is 
somewhat limited. The first two of these ratios (0/S and 
S / Al) involve elements observed in neon nova shells (Ta- 
ble [TJ, while the other two ratios (0/P and P/Al) involve 
phosphorous, which has not been observed. Thus, our re- 
sults call for improved observations and spectral analysis 
of neon nova ejecta (i.e., for phosphorous). These ratios 
are also uncertain by factors of 3-6, so future laboratory 
measurements of nuclear reactions would clearly improve 
their utility as nova thermometers. However, these four 
ratios also demonstrate the strongest temperature depen- 
dence of the potential thermometers, which mitigates the 
effect of reaction rate uncertainties. Therefore, the two 
currently observed ratios (0/S and S/Al), despite their 
large uncertainties, can be effectively used to narrow the 
range of white dwarf mass for a nova, as will be shown 
in Section El for V838 Her (Figure [6]). 

Detailed information on our results is provided in Ta- 
bleU The first four entries (N/0, N/Al, 0/S, and S/Al) 
relate to nova thermometers involving elements observed 
in neon nova shells, while the last four involve sodium 
and phosphorous which have not been observed. Col- 
umn 2 lists the maximum range of a given elemental 
abundance ratio versus peak temperature from our post- 
processing analysis. Recall that the larger this value (i.e., 
the steeper the dependence on temperature), the more 
useful a ratio will be as a nova thermometer. 

The next columns list the two reactions whose current 
rate uncertainties most significantly infiuence a given el- 
emental abundance ratio. A number of important con- 
clusions can be drawn from this information. First, the 
"second-most important reactions" (columns 5 and 6) 
give rise to very small elemental abundance ratio vari- 
ations (less than 30%), a number that is small com- 
pared to the uncertainty in observed elemental abun- 
dances (Table [1]). Second, the nova thermometers 0/S, 
S/Al, 0/P, and P/Al, which have the steepest depen- 
dence on peak temperature (by factors of 220-540; see 
column 2) also show the largest sensitivity to current 
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Figure 4. (Color online) Eight Ratios of elemental abundances (mass fractions), prime candidates for nova thermometers, that show a 
steep and monotonic dependence on peak temperature. Solid black and red lines correspond to element ratios obtained from our post- 
processing and hydrodynamic simulations, respectively Dashed lines indicate the uncertainty bands obtained by independently varying 
all relevant nuclear reaction rates within their uncertainties. Broad uncertainty bands call for future laboratory measurements of nuclear 
reactions (Table [4|). Left Panel: All ratios shown in the left column involve elements that have been observed in neon nova shells (see 
Table [T}. Right Panel: All ratios shown in the right column involve elements (Na and P) that have not been observed in neon nova shells 
(see Table [U. 
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reaction rate uncertainties (by factors of 3-6; see col- 
umn 4). Third, the latter four nova thermometers 
are mainly sensitive to current rate uncertainties of a 
single reaction, ^^P(p,7)^^S. This reaction involves a 
short-lived nuclide (t 1/2 =2.498 min) and has not been 
measured directly yet because a sufficiently intense ^^P 
beam is currently lacking. Indirect nuclear structure 
studies, in ord er to improve the r eaction rate, have 
been reported (iDohertv et al.l I2Q12I : iParikh et aLl I2Q11I : 
IWrede et al.ll2QQ9h . At present, the spin and parity as- 
signments for some of the important threshold states in 
the ^^S compound nucleus are ambiguous. In addition, 
none of the proton partial widths for these levels are 
known experimentally. Therefore, we did not estimate 
the rates using the Monte Carlo procedure described in 
Sec. 3. In the absence of more reliable information, 
we adopted t he Hauser-Feshbach (statistical model) esti- 
mate of Ra uscher fc ThielemannI (j^OOO), assuming a fac- 
tor of 10 uncertainty in the classical nova temperature 
range. Our adopted rates for this reaction agree within 
their uncertainties with those of iParikh et al.l (|2Qlll ). 
which were obtained using a different procedure. 

The influence of the ^^P(p,7)^^S rate on the abun- 
dances in the mass region above Si has already been 
pointed out bv iJose. Coc. fc Hernanzl (|2QQ1) and, more 
recently, by IParikh et al.l (|2Q11[ ). We have quantified 
its impact on nova thermometers in the present work. 
Therefore, our results provide additional motivation for 
new laboratory measurements at radioactive ion beam 
facilities. 

In general, two different techniques can be used to 
study the impact of reaction-rate variations on nucle- 
osynthesis: variation of individual rates, as done in the 
present work, and simultaneous variation of all rates us- 
ing a Monte Carlo approach. The advantage of the first 
method is that the impact of a specific reaction can be 
quantified in a straightforward manner, whereas the sec- 
ond technique allows for the influence of correlations be- 
tween interactions in the network to be taken into ac- 
count automatically. While our results indicate that pri- 
mary sources of uncertainty, such as uncertainty in the 
^^P(p,7)^^S rate, overwhelm the effects of any secondary 
source of uncertainty, we chose to verify our individual 
variation method. A recently implemented Monte Carlo 
reaction network, which included the instant-mixing and 
no- mixing approximations, was used to reproduce the 
nuclear thermometers and their uncertainties. This ap- 
proach supports the results of individual reaction rate 
variation when correla tion effects are take n into account, 
similar to the work of IParikh et al.l (j2QQ8D . 
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5. COMPARISON WITH OBSERVATION 

We will now compare the abundances from our sim- 
ulations with observations (Table [1]). A first general 
impression can be obtained from the top panel of Fig- 
ure [5l showing a three-element abundance plot involving 
N, O and Al, i.e., N/0 vs. N/Al. Our model predic- 
tions are shown as thick and thin solid, black lines (see 
Figure 2]). As already discussed in Sec. 4, the simulated 
N/0 and N/Al abundance ratios have relatively small 
uncertainties, which is reflected by the narrow region be- 
tween the two thin solid lines. These abundance ratios 
represent useful thermometers since they are not signif- 
icantly affected by current reaction rate uncertainties. 
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Figure 5. (Color online) Three-element plots used to compare el- 
emental ratios predicted by post-processing calculations to that of 
classical nova observations. The thick, solid black line corresponds 
to the elemental ratios predicted by post-processing calculations 
using recommended reaction rates. The thinner black lines de- 
note the range of uncertainty in these ratios due to reaction rate 
variations (see Figure [4|). The data points give the elemental ra- 
tios as observed in several neon novae (see Table Top Panel: 
Three-element plot of N, O, and Al. The peak temperature in the 
simulation results increases from left to right. Bottom Panel: 
Three-element plot of O, Al, and S. Note that peak temperature 
in the simulation results increases from right to left in this figure. 



The elements N, O, and Al have been simultaneously 
observed in five neon novae, according to Table [TJ and 
their abundance ratios are displayed as data points. It is 
apparent that in this case the uncertainties in the obser- 
vational data, with the exception of V382 Vel (green), far 
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Table 4 

Uncertainty Sources of Predicted Elemental Abundance Ratios in Neon Nova Shells 



Ratio 


Range ■'^ 


Primary Source 


Secondary Source 






Reaction 


Uncertainty-^ 


Reaction 


Uncertainty-^ 


N/O 


13.4 


160(p,7)17F 


1.16 


13n(p,7)14o 


1.06 


N/Al 


5.59 


20Ne(p,7)2iNa 


1.29 


13n(p,7)14o 


1.18 


O/S 


332 


30p(p,7)3iS 


3.36 


28Si(p,7)^^P 


1.09 


S/Al 


529 


30p(p,^)31S 


4.62 


28Si(p,7)29p 


1.12 


O/Na 


10.8 




1.16 


^^Ne(p,7)^-^Na 


1.12 


Na/Al 


6.83 


23Na(p,7)24Mg 


1.19 


20Ne(p,7)2iNa 


1.10 


O/P 


541 


30p(p,^)31S 


6.44 


i60(p,7)^"F 


1.26 


P/Al 


216 


30p(p,^)31S 


6.53 


20Ne(p,7)2iNa 


1.22 



^ Factor variation of final elemental abundance ratio over range of nova models 
explored in present work, obtained from the solid black lines shown in Figured 
^ Factor uncertainty of final elemental abundance (mass fraction) ratio, caused by 
varying rate of individual reaction within its current uncertainty. 



exceed those in the simulations. We conclude that the 
N/O and N/Al abundance ratios will become more useful 
thermometers for a number of novae when more reliable 
abundances can be determined from UV, IR, and optical 
spectra. A three-element abundance plot involving O, 
Al, and S, i.e., O/S vs. S/Al, is shown in bottom panel 
of Figure [5l We previously mentioned in Sec. 4 that the 
predicted O/S and S/Al abundance ratios have rather 
large uncertainties, which is reflected by the broad region 
between the two thin, solid lines. These uncertainties 
originate from the poorly known ^^P(p,7)^^S reaction 
rate (Figure |4] and Table B]). Only for a single neon nova, 
V838 Her (blue), have all three elements been observed 
simultaneously (Table [T]). It is apparent in this case that 
the observational uncertainties are relatively small com- 
pared to the stellar model results. We conclude that the 
O/S and S/Al abundance ratios will become more useful 
thermometers for neon novae when the ^^P(p,7)^^S reac- 
tion rate can be determined experimentally and Al and 
S are simultaneously observed in additional novae. 

A more detailed comparison between prediction and 
observation can be made by comparing observed abun- 
dance ratios to two or more corresponding nuclear ther- 
mometers. We are interested to find out if, for a given 
neon nova, the derived peak temperatures from two or 
more nuclear thermometers are in mutual agreement. 
Only Nova V838 Her exhibits observed elemental abun- 
dance ratios corresponding to all four currently observ- 
able nuclear thermometers (N/O, N/Al, O/S, and S/Al), 
while observed abundance ratios corresponding to two 
nuclear thermometers (N/O and N/Al or N/O and O/S) 
are available for five other neon novae (Table [T]). Thus 
we are disregarding in the following discussion neon no- 
vae V4160 Sgr and V1974 Cyg, for which only a single 
nuclear thermometer (N/O) is available. 

Nova V838 Her: Predicted and observed values for 
Nova V838 Her are presented in Figure [6l The observed 
N/O and O/S abundance ratios strongly restrict the peak 
temperature to a range of Tpeak = 0.30 — 0.31 GK, cor- 
responding to a white dwarf mass of Mwd = 1.34—1.35 
M0. For this temperature range, the observed N/Al 
abundance ratio barely misses the predicted values, al- 
though it is clear that the uncertainty of the observed 
value is too large for the observation to be of any use. 
An improved estimate of the white dwarf mass range 



could be obtained by reducing the nuclear uncertainties 
in the O/S and S/Al thermometers. 

Our above estima te supports the suggestion of 
iPolitano et al.l (|1995[ ) that the white dwarf in this 
system "is very massive." On the other hand, 
our estimate significantly exceeds the result of 
IWanajo. Ha shimot o, fc Nomotol () 19991 ). which suggests 
a value of Mwd = 1-05 Mq. Notice that our simulations 
reproduce the observed abundance ratios in the narrow 
peak temperature range shown in Figure [6l although no 
"breakout" of matter from the CNO mass range to heav- 
ier masses has occurred in any of our nova models. Thus 
we find no compelling evidence f or a breakout in V8 38 
Her, contrary to previous results (|Schwarz et al.ll2QQ7[ ). 

Nova V382 Vel: The N/O and N/Al thermometers, 
which exhibit relatively small nuclear uncertainties, re- 
strict the peak temperature to a narrow range of Tpea/c = 
0.23 — 0.24 GK, corresponding to a white dwarf mass of 
Mwd = 1-18 — 1.21 M©. For this relatively low peak 
temperature range, only a modest overproduction (fac- 
tor 2) of sulfur is expected (Figure Ej), consistent with its 
lack of observation in this neon nova. 

Nova V693 CrA: The observed N/O and N/Al abun- 
dance ratios have relatively large uncertainties and, 
consequently, the peak temperature is not well con- 
strained. The best estimate, based on these two nu- 
clear thermometers, results in an upper limit of T^eafc < 
0.28 GK, corresponding to a white dwarf mass of 
Mwd ^ 1-3 M©. An improved estimate could be 
obtained with more rehable observed abundance ra- 
tios. Our result is sig:nificant ly larger than that of 
IWanajo. Hashimoto, fc Nomotol (^1999.), who suggest a 
value of Mwd = 1.05 M©. 

Nova LMC 1990#1: The uncertainties of the observed 
N/O and N/Al elemental abundance ratios are too large 
to narrowly constrain the peak temperature range, and 
thus we only obtain an upper limit of Tp^ak < 0.24 GK, 
corresponding to a white dwarf mass of Mwd < 1.2 M©. 
Improved observations would certainly result in a more 
reliable prediction. 

Nova QU Vul: The observed N/O ratio hints at a 
low peak temperature. However, the observed N/Al 
abundance ratio is smaller than the predicted values 
over the entire peak temperature range explored in the 
present work. Therefore, we can only conclude that 
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Figure 6. (Color online) Detailed comparison between prediction 
and observation of Nova V838 Her for all four currently observable 
nova thermometers (N/O, N/Al, O/S, and S/Al). The solid black 
lines correspond to elemental ratios obtained from hydrodynamic 
simulations. The dashed lines indicate the uncertainty bands which 
were obtained by post-processing and extrapolated to the hydro- 
dynamic results. The horizontal red bands correspond to the range 
of observed values of V838 Her for each elemental ratio. The verti- 
cal red bands show the estimated peak temperature range for this 
nova based on the intersection of predicted and observed values. 



the peak temperature is T^ea/c < 0-22 GK, correspond- 
ing to a white dwarf mass of Mwd < 1-2 M©. This 
value is near the low er limit for single ONe white dwarfs 
(jDohertv et al."2010'). It also agrees with the suggestion 
of .Wanajo. Hashimoto, fc Nomoto. (1999.) . who report a 
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Figure 7. (Color online) Peak temperature and white dwarf mass 
ranges as given by intersection of observed and predicted values of 
nova thermometers (see, for example, Figure O. Bracketed bars 
refer to white dwarf mass ranges where a definite range could be 
established. Bars terminating with an arrow refer to white dwarf 
mass ranges where only an upper limit could be determined. 

range of Mwd = 1-05 - 1.1 M©. 

Nova VI 065 Cen: In this case, observed N/O and O/S 
abundance ratios are available. However, no range of 
peak temperatures can be found that give consistent re- 
sults for these two nuclear thermometers. 

The resulting peak temperature and white dwarf mass 
ranges for V838 Her, V382 Vel, V693 CrA, LMC 1990#1, 
and QU Vul are shown in Figure [71 Although these re- 
sults are encouraging, it must be kept in mind that the 
present simulations are obtained for one particular choice 
of accretion rate, white dwarf luminosity, and mixing 
fraction between accreted and white dwarf matter. Evo- 
lutionary sequences for other choices of these parameters 
are in progress and will be presented in a forthcoming 
pubhcation. 

6. SUMMARY AND CONCLUSIONS 

One goal of the present work was to investigate how 
useful elemental abundances are for determining the peak 
temperature, which is strongly correlated with the un- 
derlying white dwarf mass, achieved during neon nova 
outbursts. Another aim was to determine how robust 
such "nova thermometers" are with respect to currently 
uncertain nuclear physics input. To this end, we first 
presented updated observed abundances for several neon 
novae, and we performed new hydrodynamic simulations 
of neon novae, with peak temperatures in the range of 
228 MK to 313 MK. We found that the most useful ther- 
mometers (i.e., those elemental abundance ratios with 
the steepest monotonic dependence on peak tempera- 
ture) are N/O, N/Al, O/S, S/Al, 0/Na, Na/Al, 0/P, 
and P/Al. These vary by factors of 13, 6, 330, 530, 11, 
7, 540, and 220, respectively, over the temperature range 
explored here (see Table |4]). 

Next, we investigated the sensitivity of these nova ther- 
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mo meters to thermonuclear reaction rate variations. The 
ratios N/0, N/Al, 0/Na, and Na/Al are robust, in the 
sense that they are affected by uncertain reaction rates 
by less than 30%. At present, only the ratios N/0 and 
N/Al are useful thermometers since these elements have 
been observed in neon novae (see Table [T]). The last 
two ratios, 0/Na and Na/Al, are not currently useful as 
sodium has not yet been detected in classical nova spec- 
tra. We also find that for these four elemental abundance 
ratios the dependence on peak temperature (factor vari- 
ations of 6-13; Table llj is not as strong as for some other 
ratios, mentioned below. 

The thermometers 0/S, S/Al, 0/P, and P/Al reveal 
far steeper monotonic dependences on peak temperature, 
by factors of 220-540. However, their current drawback 
is the strong dependence on the uncertain ^^P(p,7)^^S 
reaction rate. Thus, our study provides additional mo- 
tivation for new laboratory measurements of this crucial 
nuclear reaction. In addition, the last two of these ra- 
tios, 0/P and P/Al, involve phosphorus, another element 
for which reliable abundances have not yet been derived 
from classical nova spectra in outburst. 

Finally, we compared our model predictions to elemen- 
tal abundances observed in neon nova shells. Based on 
the present nuclear thermometers, we obtain the follow- 
ing estimates for the underlying ONe white dwarf masses: 
1.34-1.35 Mo (V838 Her), 1.18-1.21 M© (V382 Vel), <1.3 
Mo (V693 CrA), <1.2 Mo (LMC 1990#1), and <1.2 Mo 
(QU Vul). Presently no white dwarf mass range based 
on elemental abundances can be derived for nova VI 065 
Gen. These predictions could be improved if more reli- 
able abundances from the analysis of UV, IR, and opti- 
cal spectra become available in the future. The present 
simulations were obtained for one particular choice of ac- 
cretion rate, white dwarf luminosity, and mixing fraction 
between accreted and white dwarf matter. Evolutionary 
sequences for other choices of these parameters will be 
presented in a forthcoming publication. 

We would like to thank Anne Sallaska, Greg Schwarz, 
Steve Shore, and Karen Vanlandingham for fruitful dis- 
cussions. This work was supported in part by the US De- 
partment of Energy under grant DE-FG02-97ER41041, 
the National Science Foundation under award number 
AST-1008355, the Spanish MIGINN grants AYA2010- 
15685 and EUI2009-04167, the E.U. FEDER funds, the 
ESF EUROGORES Program EuroGENESIS, and by 
Arizona State University funding from the National Sci- 
ence Foundation and the National Aeronautics and Space 
Administration. 

REFERENGES 

Althaus, L. G., Corsico, A. H., Isern, J., &; Garcia-Berro, E. 2010, 

A&A Rev., 18, 471 
Andrea, J., Drechsel, H., & Starrfield, S. 1994, A&A, 291, 869 
Bode, M., & Evans, A. 2008, Classical Novae (2nd ed.; 

Cambridge, UK: Cambridge University Press) 
Casanova, J., Jose, J., Garcia-Berro, E., Shore, S. N., Calder, A. 

C. 2011, Nature, 478, 490 
Doherty, C. L., Siess, L., Lattanzio, J. C, &; Gil-Pons, P. 2010, 

MNRAS, 401, 1453 
Doherty, D. T., Lotay, G., Woods, P. J., Seweryniak, D., 

Carpenter, M.P., Chiara, C. J., David, H. M., Janssens, R. V. 

P., Trache, L., & Zhu, S. 2012, Phys. Rev. Lett., 108, 262502 



Ferland, G.J. 2003, ARA&A, 41, 517 

Gehrz, R. D., Truran, J. W., WiUiams, R. E., & Starrfield, S. 

1998, PASP, 110, 3 
Helton, L., Woodward, C. E., Walter, F. M., et al. 2010, A J, 140, 

1347 

Hernanz, M. 2008, in: Classical Novae (2nd Ed.), Eds. M.F. Bode 

and A. Evans (Cambridge University Press) 
Iliadis, C, Champagne, A., Jose, J., Starrfield, S., &; Tupper, P. 

2002, ApJS, 142, 105 
Iliadis, C. 2007, Nuclear Physics of Stars (Weinheim: Wiley- VCH) 
Iliadis, C, Longland, R., Champagne, A. E., Coc, A., & 

Fitzgerald, R. 2010a, Nucl. Phys. A, 841, 31 
Ihadis, C, Longland, R., Champagne, A. E., & Coc, A. 2010b, 

Nucl. Phys. A, 841, 251 
Iliadis, C, Longland, R., Champagne, A. E., &; Coc, A. 2010c, 

Nucl. Phys. A, 841, 323 
Iliadis, C, Champagne, A., Chieffi, A., &; Limongi, M. 2011, 

ApJS, 193, 16 
Jose, J. & Hernanz, M. 1998, ApJ, 494, 680 
Jose, J., Coc, A., & Hernanz, M. 2001, ApJ, 560, 897 
Jose, J., Hernanz, M. & Ihadis, C. 2006, Nucl. Phys. A, 777, 550 
Jose, J., & Shore, S.N. 2008, in: Classical Novae (2nd Ed.), Eds. 

M.F. Bode and A. Evans (Cambridge University Press) 
Kovetz, A. & Prialnik, D. 1997, ApJ, 477, 356 
Fodders, K., Palme, H., &; Gail, H.-P. 2009, Landolt-Bornstein 

New Series VI/4B (Springer 2009) 
Longland, R., Iliadis, C, Champagne, A. E., Newton, J. R., 

Ugalde, C, Coc, A., &; Fitzgerald, R. 2010, Nucl. Phys. A, 841, 

1 

Mason, E. 2011, A&A, 532, Lll 

Parikh, A., Jose, J., Moreno, F., & Iliadis, C. 2008, ApJS, 178, 
110 

Parikh, A., Wimmer, K., Faestermann, T., et al. 2011, 

Phys. Rev. C, 83, 045806 
Politano, M., Starrfield, S., Truran, J. W., Weiss, A., &; Sparks, 

W. M. 1995, ApJ, 448, 807 
Prialnik, D. 1986, ApJ, 310, 222 

Rauscher, T. & Thielemann, F. K. 2000, ADNDT, 75, 1 
Ritossa, C, Garcia-Berro, E., & Iben, Jr., I. 1996, ApJ, 460, 489 
Schwarz, G. J. 2002, ApJ, 577, 940 

Schwarz, G. J., Shore, S. N., Starrfield, S., &; Vanlandingham, K. 

M. 2007, ApJ, 657, 453 
Shara, M. M., Yaron, O., Prialnik, D., Kovetz, A., &; Zurek, D. 

2010, ApJ, 725, 831 
Shore, S. N., et al. 2003, AJ, 125, 1507 
Shore, S. N. 2012, private communication 

Smith, M. S., Hix, W. R., Parete-Koon, S., Dessieux, L., Guidry, 

M. W., Bardayan, D. W., Starrfield, S., Smith, D. L., & 

Mezzacappa, A. 2002, Proceedings of the International 

Conference on Classical Nova Explosions, eds. M. Hernanz and 

J. Jose, American Institute of Physics 
Snijders, M. A. J., Batt, T. J., Roche, P. F., Seaton, M. J., 

Morton, D. C, Spoelstra, T. A. T., &; Blades, J. C. 1987, 

MNRAS, 228, 329 
Starrfield, S., Sparks, W. M., &; Truran, J. W. 1986, ApJ, 303, L5 
Starrfield, S., Truran, J. W., Wiescher, M. C, & Sparks, W. M. 

1998, MNRAS, 296, 502 
Starrfield, S., Sparks, W. M., Truran, J. W., &; Wiescher, M. C. 

2000, ApJS, 127, 485 
Starrfield, S., Ihadis, C, & Hix, W. R. 2008, in: Classical Novae 

(2nd Ed.), Eds. M.F. Bode and A. Evans (Cambridge 

University Press) 
Starrfield, S., Ihadis, C, Hix, W. R., Timmes, F. X., &; Sparks, 

W. M. 2009, ApJ, 692, 1532 
Vanlandingham, K. M., Starrfield, S., & Shore, S. N. 1997, 

MNRAS, 290, 87 
Vanlandingham, K. M., Starrfield, S., Shore, S. N., & Sonneborn, 

G. 1999, MNRAS, 308, 577 
Vanlandingham, K. M., Schwarz, G. J., Shore, S. N., Starrfield, 

S., &; Wagner, R. M. 2005, ApJ, 624, 914 
Vanlandingham, K. M. 2012, private communication 
Wagner, R. M. &; DePoy, D. L. 1996, ApJ, 467, 860 
Wanajo, S., Hashimoto, M., &; Nomoto, K. 1999, ApJ, 523, 409 
Wrede, C, Caggiano, J. A., Clark, J. A., Deibel, C. M., Parikh, 

A., &; Parker, P. D. 2009, Phys. Rev. C, 79, 045803 
Yaron, O., Prialnik, D., Shara, M. M., & Kovetz, A. 2005, ApJ, 

623, 398 



